Multiparticle entanglement purification for graph states 



W. Dur, H. Aschauer and H.-J. Briegel 
Sektion Physik, Ludwig-Maximilians-Universitat Miinchen, Theresienstr. 37, D-80333 Miinchen, Germany. 

(Dated: February 1, 2008) 

We introduce a class of multiparticle entanglement purification protocols that allow us to distill a 
large class of entangled states. These include cluster states, GHZ states and various error correction 
codes all of which belong to the class of two-colorable graph states. We analyze these schemes under 
realistic conditions and observe that they are scalable, i.e. the threshold value for imperfect local 
operations does not depend on the number of parties for many of these states. When compared 
to schemes based on bipartite entanglement purification, the protocol is more efficient and the 
achievable quality of the purified states is larger. As an application we discuss an experimental 
realization of the protocol in optical lattices which allows one to purify cluster states. 

PACS numbers: 03.67.-a, 03.67.Mn, 03.67.Pp 



Entangled states of multi-partite systems are expected 
to play a significant role for applications in quantum com- 
putation and quantum communication. Quantum error- 
correcting codes, as the most prominent example, are 
highly entangled states which are used to encode quan- 
tum information and to protect it against the damaging 
effects of decoherence. It has been shown that similar en- 
tangled states may be used in multi-party communication 
scenarios beyond teleportation and quantum key distri- 
bution fjj^uch as secret sharing or secure function eval- 
uation 0, Q , but also in more practical applications such 
as the improvement of existing frequency standards 0. 
In the context of quantum computation, specific multi- 
partite entangled states — the so-called cluster states |5j 
- have even been shown to constitute a universal resource 
for measurement- based quantum computation . These 
states can be created e.g. via an Ising interaction between 
neighboring particles on a lattice. A specific realization 
of such a system is based on neutral atoms in an optical 
lattice where, starting from a Mott-insulator state 0, 
cluster states could be created by a simple interferomet- 
ric process 0, . 

When talking about entanglement, it is mandatory to 
consider the effect of decoherence. Finite temperature (in 
the case of interacting particles) or a noisy communica- 
tion channel (in the case of distributed quantum systems) 
lead to decoherence in the state space of the information 
carriers, with a rate that increases typically in propor- 
tion with the number of particles. Therefore only mixed 
states rather than pure states will be available, and the 
achievable fidelity is expected to decrease exponentially 
both with time and with the size of the system. It is thus 
not clear whether an entangled state of a large number 
of particles can be created and maintained in practice. 

Standard methods have been developed to stabilize 
quantum states against the detrimental effects of deco- 
herence and noise. Quantum error correction is one gen- 
eral method, but the acceptable noise level for the quan- 
tum operations is extremely small (of the order of 1CP 4 ) 
0. Entanglement purification [Tol lllj| is an alternative, 
but protocols exist only for the purification of a special 
type of states, namely states which are equivalent, up 



to local unitary operations, to the so-called "cat" states 
10)®* + 11)®* [la [II El 111. Furthermore, it needs to 
be investigated whether these protocols are at all appli- 
cable under realistic conditions. While in the bipartite 
case certain entanglement purification protocols turn out 
to be remarkably robust against the influence of imper- 
fect local operations - allowing for threshold values at 
the order of several percent 01 - > it is n °t clear whether 
this is also true in the multipartite case. In particular, 
it is not known how the threshold for the required preci- 
sion of local control operations scales with the size of the 
system, a crucial question if one wants to establish and 
maintain quantum entanglement in larger systems. 

In this paper, we introduce a class of entanglement pu- 
rification protocols, which can be used to purify a large 
class of multi-partite entangled states, in particular all 
two-colorable graph states. These states include, for ex- 
ample, various quantum error correcting codes [15|, the 
generalized GHZ states (or cat states), and the cluster 
states, all of which are resources for the applications men- 
tioned above. We analyze these protocols under realis- 
tic conditions, i.e. when the local operations required in 
the purification process are imperfect. We find that: (i) 
The value of the purification threshold, i.e. the acceptable 
noise level for local operations under which the protocol 
still purifies, does not depend on the number of particles, 
but is influenced by the maximal degree of the graph as- 
sociated with the state, (ii) The described multi-particle 
entanglement purification protocols (MEPP) are gener- 
ally not only more efficient than any method based on bi- 
partite entanglement purification protocols (BEPP), but 
in the case of imperfect local operations the achievable 
fidelity is even higher, (in) The entanglement of the pu- 
rified multipartite states is private |l6j. thus providing a 
secure resource for applications in multiparty communi- 
cation scenarios 3]. We close by describing how MEPP 
could be implemented experimentally in optical lattices, 
allowing one to increase the fidelity of cluster states. 

We start out by introducing two-colorable graph 
states, and purification protocols which allow one to pu- 
rify all states of this kind. Consider a graph G = (V, E) 
which is a set of vertices V connected in a specific way 
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by edges E specifying the neighborhood relation between 
vertices. With this graph we associate N = \V\ com- 
muting correlation operators Kj — <7i ]J^ k -^ eE erg . 
The joint eigenstates of Kj, 1*^^...^), p,j G {0,1}, 
the graph states associated with G, fulfill the eigenvalue 
equation Kj\^^ 2 ... m ) = (-l)^|* MlMa ...^) Vj. Note 
that {|^ lM2 ... AlJV )} G form a basis in H = ((F 2 )® N . We 
restrict our attention to two-colorable graphs, that are 
graphs for which a partition of the vertices into two dis- 
joint sets V a UVb — V exists such that no vertices within 
a set are connected by edges. The states arising from 
such two-colorable graphs, which we call two-colorable 
graph states (TCGS), include various interesting MES, 
ranging from iV-particle GHZ states over codewords for 
quantum error correction codes to cluster states. In what 
follows, we will introduce protocols which allow one to 
purify arbitrary TCGS. That is, given M copies of an 
arbitrary TV-partite mixed state p, we establish for each 
two-colorable graph G a local protocol which is capable 
of creating the pure state |^g)G as output state, pro- 
vided the initial state p fulfills certain requirements (e.g. 
has sufficiently high fidelity with respect to \^^)g)- 

Let us consider an arbitrary but fixed two-colorable 
graph G with vertices V = V A UV B , N A = \V A \,N B = 
\Vb\ and N = N A + Nb spatially distinct parties each 
holding one of the N particles belonging to a gen- 
eral mixed state p. One can depolarize the state p to 
a state pa which is diagonal in the graph-state basis 
{|^ , juijU2..-auv}}g without changing the diagonal elements 
by a probabilistically applying the local operations cor- 
responding to the operators Kj [17\. That is, without 
loss of generality we can consider mixed states diagonal 
in the graph-state basis, 

PG= zl X Pa,Pb\^Pa-Pb)(^Pa.Pb\> (i) 

Pa,Pb 

where we introduced the multi indices ft Al p B with e.g. 
pZ A = /iji/iia . . . t L i NA indicating all vertices in the set V A . 
We now introduce two protocols P1,P2 acting on two 
identical copies pi = P2 = Pg, Pyi = Pi ® P2- In protocol 
PI, all parties which belong to set V A [V B ] apply local 
CNOT operations ^3 to their particles, with the particle 
belonging to p\ \pi\ as source, pi [pi] as target. From the 
eigenvalue equations for \^p A .p B ) follows that the action 
of the local CNOT operations is given by 

\^Pa,Pb ) V^va.vb ) > \^ Pa,Pb(BPb ) \^ va®Pa,vb )i (2) 

where fi A (Bv A denotes bitwise addition modulo 2. For in- 
stance, if jl A = Pip.3^5, fi A ®V A = fll®Vl, fl3®V3, fl5®V5. 

All particles of pi belonging to set V A [Vb] are then 
measured in the eigenbasis of a x [a z ] respectively, yield- 
ing results (-1)^ [(-l) Cfc ] with £j,( k G {0,1}. The 
first state is only kept if the measurement outcomes 
fulfill (£j + J2{k j}£E Cfe) m °d2 = Vj which implies 
p A © v A — 0. Eq. @ makes it evident that in this proce- 
dure information about the first state p\ - encoded into 



fi A ,pB ~ is transferred to the second state p2 and re- 
vealed by the measurement. This is in analogy to the 
standard recurrence protocols [13, m ei m and it is 
the key point to purifying mixed states. 

The resulting state p after this procedure is again 
diagonal in the graph-state basis, with new coeffi- 
cients XjA-,jB — H{{u b ,Pb)\vb®Pb=1b} 2K^1a,v b ^ia,Pbi 

where if is a normalization constant such that tr(/5) = 
1. In the protocol P2, the role of sets V A 
and Vr is interchanged, which leads to A'- - = 

T,{{u a ,Pa)\Va®Pa=ia} m X tA,iB^p A ,iB- The total purifi- 
cation protocol corresponds to an iterative application 
of subprotocols PI and P2, always using two identical 
copies of the multiparticle states obtained in the previ- 
ous round as input states. For GHZ states, this protocol 
is equivalent to the one introduced in Ref. [T^ . and fur- 
ther analyzed in Ref. [13(, and may thus be viewed as a 
generalization to the class of all TCGS. 

To gain analytical insight into this procedure, we con- 
sider the example of rank 2 Na mixed states of the form 
PA = Ep A ^Wl^/U.ffX^.Bl- The application of pro- 
tocol PI leads to another state of this form with new 
coefficients A,- s = \1 -.IK and K = V- \ 2 _ 

MA,0 Pa,0' ^HA p A ,0 

That is, the largest coefficient is amplified with respect 
to the other ones and iteration of the protocol allows 
one to produce pure graph states. Consider for in- 
stance the one parameter family p A (F) with Agg = F, 

A- A g = (1 - F)/{2 Na - 1) for pi A ^ 0, where F is 
the fidelity of the desired state. Application of PI pre- 
serves the structure of those states and leads to F = 
F 2 /[F 2 + (1 - F) 2 /{2 Na - 1)]. This map has F = 1 as 
attracting fixed point for F > 1/2 Na . In these exam- 
ples, application of protocol PI alone is sufficient as only 
information about p, A has to be extracted. 

For general full rank mixed states, both PI and P2 will 
have to be applied in order to reveal information about 
Pa [Pb] respectively. While PI increases the weight of 
coefficients Ag - b , P2 amplifies coefficients A~ g which 
together leads to the amplification of Ag g given the ini- 
tial fidelity is sufficiently high. The action of PI, P2 
correspond to non-linear mappings of a large number of 
independent variables (in total 2 N — 1), which makes an 
analytic treatment of the purification protocol very diffi- 
cult. We have however analyzed the protocol numerically 
and found that it is indeed capable to purify noisy TCGS 
arising from various kinds of noise. In the following, we 
will concentrate on GHZ states and ID cluster states. 
The graph corresponding to the TV-particle GHZ state 
is given by V A = {1},V B = {2,3,...,iV} with edges 
{l,k},k G Vb, while ID cluster states correspond to a 
graph where vertex k is connected to vertex k + 1 , Vfc 
and sets V A [Vb] are formed by all odd [even] vertices 
respectively. Our results are however not restricted to 
these specific graphs. 

We consider noisy TCGS arising naturally in a mul- 
tipartite scenario where each of the N particles consti- 
tuting j^g) is subject to decoherence, arising e.g. from 
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sending the particles through noisy quantum channels. 
We consider depolarizing channels with noise parameter 
q described by 



£ k p = qp+(l- q)/21 k <g> tr k (p), 



(3) 



where in this case the channel is acting on particle k. 
The resulting multipartite state is of the form p(q) = 
£\&2 ■ ■ . £^jv l^o) (^o I wri ich we take as input states for 
our MEPP. It turns out that the tolerable amount of 
white noise per particle such that our MEPP can still 
be successfully applied is for linear cluster states essen- 
tially independent of the number N of particles, while 
for GHZ states this value decreases exponentially with 
N (see Fig. QJi)- This implies that, for the GHZ state, 
the required quality of the channel over which particles 
are distributed, such that the resulting state is still purifi- 
able, increases with N. Note that this different behavior 
is not reflected in the minimum fidelity F m - ln which in 
both cases decreases exponentially (Fig.^Ji). We remark 
that for all states p(q) where q is larger than the threshold 
value, the protocol successfully converges towards the de- 
sired state l^g). For more general states, conditions for 
convergence are difficult to determine due to the large 
number of parameters, even though convergence itself 
can straightforwardly be checked numerically. Also for 
noisy states of the form p(x) — x\*S> q) (9 $\ + (1 — x)/2 N l, 
i.e. mixtures of the desired state with a completely depo- 
larized state, the situation is similar, i.e. P m i n decreases 
exponentially with N. 
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FIG. 1: (a): Minimal value of fidelity F m m + [o] and pa- 
rameter g m j n x [*] for linear cluster states [GHZ states] for 
different number of particles A^ and perfect local operations, 
(b): Achievable fidelity of a linear cluster state with A" = 4 
using direct MEPP (solid line) and conservative upper bound 
for methods based on BEPP (dashed line) for different errors 
in local operations p. 

Until now we have assumed that local operations - 
in particular CNOT operations - are perfect. In prac- 
tice, however, these operations as well as measurements 
will be imperfect. We now investigate the influence of 
errors in the local operations on the MEPP. We will 
consider a simple error model where local two-qubit op- 
erations are described by the completely positive map 
£u jk P = Ujk[£j£kp]Uj k , where £ k ,£j are given by Eq. © 
with error parameter p. That is, an imperfect operation 
is described by first applying local white noise with prob- 
ability (1 — p) independently on the qubits, followed by 



the perfect unitary operation. We have also investigated 
more general error models, e.g. correlated white noise, 
and errors in the measurement process, observing essen- 
tially the same behavior as for this simple model. 

We have investigated the dependence of the minimal 
required fidelity and the maximal reachable fidelity for 
linear cluster states of different length on the error pa- 
rameter p (see Fig. as well as the threshold value 
p m i n until which our MEPP can be successfully applied 
(Fig 03 ) . As can be seen from Fig. Pmin is almost 
independent of the number of particles N, it even seems 
that for larger A?" the tolerable amount of noise per oper- 
ation is larger. For GHZ states, in contrast, the thresh- 
old value Pmi n increases with N, i.e. the requirements 
to purify GHZ states with larger N are more stringent. 
The qualitatively different behavior of p m in for these two 
TCGS can be understood within a restricted error model. 
We consider noisy operations where only particles in Vb 
are subjected to bit-flip errors - described by the map 
£kP = VP + (1 — P)/2{p + o- x pa x ) - and entanglement 
purification is performed by applying protocol PI. For 

GHZ states we have that equivalent to phase 

flip errors in particle 1, cr^ 1 Vfcs E Vb as can be seen 
from the eigenvalue equations, while for closed linear 
cluster states a x 2k ^ corresponds to eri 2 ' c ~ 1 ^(7i 2fc+1 \ Thus 
for GHZ states, errors at N — I particles affect particle 
I and the errors accumulate. For linear cluster states, 
however, each particle in Va is only affected by errors at 
two other particles in Vb, independent of N. One can 
show analytically that even for this restricted class of er- 
rors, for GHZ states p m in increases exponentially with 
N, p m in = 1/2 1 '( N ~ 1 '. For linear cluster states one finds 
that for input states of the form pa{F) and a fixed er- 
ror rate, e.g. p — 0.8, one application of the protocol PI 
increase the fidelity F = x + (1 — x)/2 N ^ 2 in the range 

2 -0.337V < x < 2 -0.009JV [f^ That for each N there 

exists a finite regime where entanglement purification is 
possible and the threshold value p m in converges for large 
N towards p min » 0.4938. 
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FIG. 2: (a): Maximal reachable fidelity F ma x and minimal 
required fidelity P m i n plotted against error parameter p (lo- 
cal operations) for density operators arising from single-qubit 
white noise. Curves from top to bottom correspond to linear 
cluster states with N — 2, 4, 6 particles, (b): Threshold value 
for errors in local operations p m i n for GHZ states (*) and lin- 
ear cluster states (x) with different number of particles N. 

We have also performed numerical simulations for var- 
ious other two-colorable graphs corresponding e.g. to 2D 
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cluster states and simple CSS codes [15J , observing a sim- 
ilar behavior as for linear cluster states. Since graphs cor- 
responding to a concatenation of such CSS codes are also 
two-colorable |17|. our approach may be used to purify 
entire encoding circuits used in fault tolerant quantum 
computation. While in the general case the structure of 
the graph - in particular its maximum degree - influences 
the dependence of the threshold value p m i n of the corre- 
sponding purification protocol, the exact dependence of 
Pmin on the graph is still an open problem. 

Let us compare this MEPP with schemes based on 
bipartite entanglement purification. By this we mean 
a scheme which produces bipartite entangled states by 
means of BEPP, which are then used to create a MES 
by some means, e.g. by teleportation. Even for the best 
known BEPP with respect to maximal reachable fidelity 
under imperfect local operations - which is the proto- 
col introduced in Ref. [Ill ] - and under the conservative 
assumption that the local operations involved in the cre- 
ation of multiparticle entangled states from (noisy) en- 
tangled pairs are itself error free, the maximal reachable 
fidelity F^^. for our MEPP is considerably larger than 
the upper bound for methods based on BEPP, as can be 
seen m Fig. H1d. Note that the upper bound in the lat- 
ter case solely depends on the fixed point of the BEPP. 
This shows that there are tasks that can be solved with 
multi-lateral purification (i.e. the efficiency of the pro- 
cess is non-zero) but not with bi-lateral purification (i.e. 
the efficiency is zero), where the efficiency is given by the 
ratio of surviving copies with the desired fidelity to the 
number of initial copies (see Ref. 01 f° r details). 

Finally, we propose an experimental realization of 
MEPP with neutral atoms in optical lattices. We show 
that MEPP can be used in such systems to increase the 
fidelity of cluster states. In particular, we consider the 
purification of ID cluster states in a 2D lattice. Remark- 
ably, even if the same imperfect operations are involved 
in the creation of the cluster state and in the purifica- 
tion process, MEPP allows one to enhance the achiev- 
able fidelity of the state ^?|]- For the experiment, con- 



sider a two-dimensional N x N optical lattice filled with 
one atom per lattice site. Internal states of the atoms 
- which constitute the qubits - can be manipulated by 
means of laser pulses. Interactions between neighboring 
atoms take place e.g. by state-selectively shifting the lat- 
tice, leading to a state dependent collisional phase aris- 
ing from controlled cold collisions The interaction 
Hamiltonian describing a lattice shift in the ir-direction 
is given by H x = g(t) E M (l + ^ M) )/2® (1 ~4 k+1 ' l) )/2, 
where (k,l) labels the (x, y)-coordinate of the atom and 
g(t) is a time dependent coupling parameter. Note that 
for J g(t)dt — 7r, such an interaction produces N copies of 
one dimensional cluster states along the x-direction of the 
lattice when applied to states of the form (|0) + 11))®^. 
These states can then be purified by using lattice shifts 
along the y-direction: By applying local laser pulses 
to the individual atoms, the resulting interaction for 
J g(t)dt — 7r can be converted into the proper CNOT 
operations between pairs of atoms, thereby realizing J5J). 
Using a total of two lattice shifts, this process can even 
be parallelized such that N/2 pairs of cluster states are 
purified simultaneously. The resulting states after a suc- 
cessful measurement can further be purified by applying 
protocol P2 in a similar way. 

In this paper we have introduced entanglement purifi- 
cation protocols which are capable of purifying a wide 
range of multiparticle states, namely all TCGS. In the 
case of noisy operations, these multiparticle entangle- 
ment purification protocols are not only remarkably ro- 
bust, but they allow one to reach higher fidelities than 
known schemes based on bipartite purification. We also 
found that for cluster states the threshold value for im- 
perfect local operations is independent of the number of 
particles. Our results are a step towards practical appli- 
cations based on multiparticle entangled states. 
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